Single quantum realization of a collision of two Bose— Einstein condensates 
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We propose a method for simulating a single realization of a collision of two Bose-Einstein con- 
densates. Recently in Zin et al. (Phys. Rev. Lett. 94, 200401 (2005)) we introduced a quantum 
model of an incoherent elastic scattering in a collision of two counter-propagating atomic Gaussian 
wavepackets. Here show that this model is capable of generating the data that can be interpreted 
as results of a single collision event. We find a range of parameters, including relative velocity, 
population and the size of colliding condensates, where the structure of hallo of scattered atoms in 
a single realization strongly differs from that averaged over many realizations. 
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According to the standard interpretation of quantum 
mechanics the single particle wavefunction ip(r,t) de- 
scribes probability ampUtude and ['(/'(r, i)p - probabil- 
ity density for finding this particle around r at time t. 
To make a comparison between theory and experiment 
one needs to repeat the same measurement many times 
since the agreement can be checked only on the statis- 
tical level. Recent developments in the physics of Bose 
Einstein condensates of neutral atom gasses created 
an extremely interesting situation where already a single 
measurement, due to large number of particles involved, 
reveals some statistical aspects. One prominent example 
is given by two overlapping independent Bose-Einstein 
condensates (each of them in the Fock state), where ev- 
ery single measurement of the system reveals an inter- 
ference pattern. At the same time however, no pattern 
is present in the average over many measurements 0, 1^ ■ 
This example shows that even for many particle systems 
there is a substantial difference between results of a single 
experiment and the statistical average (over many repe- 
titions of the experiment). The latter, in analogy with 
the single particle case, is much easier to handle theoreti- 
cally, the former, more straightforward for experimental- 
ist, is a real challenge for theoretician. Nevertheless, in 
some cases theoretical solutions are possible. Javanainen 
and Yoo 0| were able to prove that indeed interference 
fringes are present in single measurement of overlapping 
independent condensates. Bach and Rzazewski stud- 
ied quantum correlations in atomic systems, and noted 
that single measurements data contain information about 
these correlations. Dziarmaga proved that any Bo- 
goliubov vacuum can be brought to a diagonal form in 
a time-dependent ortho normal basis. This diagonal form 
is tailored for simulations of quantum measurements on 
an excited condensate. They illustrated their method 
using example of the phase imprinting of a dark soliton. 
Here we examine a similar approach for the case when 
the substantial depletion of the condensate occurs. 

Bragg diffraction divides a Bose-Einstein condensate 
into two initially overlapping components, which imme- 
diately start moving away from each other with high 



relative momentum. Elastic collisions between pairs 
of atoms from distinct wave packets cause losses that 
can significantly deplete the condensate. Recently |^ 
01 we introduced the quantum model of two counter- 
propagating atomic Gaussian wave packets to study such 
collisions and the transition from spontaneous to stimu- 
lated regime. Within this model correlation functions of 
all orders are accessible. It is then not surprising that 
it allows for a next conceptual step - reconstruction of 
the result in a single realization of the experiment. The 
proposed method is valid in the regime of bosonic stimu- 
lation and we show the relation to the previous analytical 
methods 0, but also how it can be used to test the va- 
lidity of stochastic method 0. 

We consider a half-collision of two Bosc-Einstcin con- 
densates initiated by Bragg scattering. Hamiltonian of 
the system in second quantization reads, see j^l2| 

H - / dV^'t(r,t)— — *(r,t) 
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y" dV^'t(r,t)#t(r^t)^'(r,i)*(r,t), (1) 



where ^'(r, i) is the bosonic field operator. Atoms in- 
teract via two-body contact interaction, determined by 
the single coupling constant g. As the system con- 
sists of two counter-propagating, highly occupied atomic 
wave-packets and the "sea" of unoccupied modes, in the 
spirit of Bogoliubov approximation we decompose the 
field operator into c— number condensate wave-function 
ipQ{r,t) -f ip^Q{r,t) and the quantum field of scattered 
atoms, (5(r, t), : 



*(r,t) = VQ(r,0+V'-Q(r,t)+<5(r,i). 



(2) 



The subscripts ±Q denote the mean momentum of the 
colliding wavepackets. In what follows we assume that 
the condensate wave-function satisfies independently the 
time-dependent GPE. We insert Eq. Q into the Hamil- 
tonian ^ and keep only terms, that lead to creation or 
annihilation of a pair of particles in the S field of scattered 
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atoms, 



H = - / d'rSHr,t)%^6iv,t) (3) 
J 2m 

+ g f dV5t(r,i)^t(r,t)V'Q(r,i)V-Q(r,i)+H.c. 



Further simplification is possible in the regime when the 
collisional time tc = a/{hQ/m) is much shorter than 
the linear disspersion time ^ld = ma"^ /h and the nonlin- 
ear disspersion time tjvn = \/ ir^^'^ma^ / gN Here a 
and iV/2 is the radius and number of particles in each of 
the colliding partners. In this regime neither dispersion 
nor the nonlinearity produce any noticeable effects on 
the condensate part during the collision and the change 
of the condensates' shape can be neglected. To model 
the condensates we use spherically symmetric Gaussian 
wave- functions 
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where r = {xi,X2,X3). In this approximation the 
Heisenberg equation for the dimensionless field opera- 
tor 6{r,t) a^/^ exp {i(3t/2) <5(r,f) in the dimensionless 
units (t/tc t, Xija — > xi (for i — 1,2,3), and with 
parameters tLo/tc = (3 and {Ild /iND)"^ — a) reads 



z/5at<5(r,<) = -i(v2+/32)5(r,t)+, 



(5) 

We decompose the field operator into a basis modes of 
spherically symmetric harmonic oscillator 9] , 

(t), (6) 

n,l,m 

as discussed in , where we solved the system of linear 
equations for operators &n,z.m(i). 

Here, following the ideas introduced in we take 
the advantage of the fact that the Hamiltonian of the 
system is quadratic in creation and anihilation operators 
and reduce the problem to independently evolving modes 
that experience squeezing. To do so we map all three 
quantum numbers n,l,m into single index ([n, / , m] — > [i]) 
and write the solution of the evolution equation @ in a 
compact form 



b{t) = Cb{0) + SP{0), 



(7) 



where vector b = (61,62,...)"^ and C and S are time 
dependent (evolution) matrices. According to [l3|, C and 
S can be written as C = UCdV' and S — U SdV^, were 
U and V are unitary and Cd and Sd are diagonal with 
non-negative eigenvalues. Hence it is possible to define 



new initial modes (vectors) a(0) = V^b{0) and new final 
modes: a{t) = Wb{t) that are related by simple diagonal 
transformation 



d{t) ^CDa{Q) + Soa^Q) 



(8) 



— sj — 1. The choice 



with a condition — 
of the time t is crucial as it determines both initial and 
final sets of modes in Eq. (jHJ. It follows from Eq. JS)) that 
the system is in the multimode squeezed state J^Oj, with 
atomic field operator 



(9) 



where ai{t) = Ciai(O) -I- SiaJ(O). Notice that the density 
matrix is these modes diagonal |0 

p(ri,r2,t) = ((5t(ri,t),5(r2,t)) =^V;(ri,t)V'.(r2,i)s-. 

i 

To find the set of -ipi (r, t) we first solve the linear equation 
© , obtain the general solution of the form Q and eval- 
uate the density matrix. Finally, following reference |3, 
we find the squeezed modes by diagonalization of the 
density matrix and the anomalous matrix. 

For highly populated modes [si ^ 1) we have Ci ~ Si. 
Then the anihilation operators can be rewritten as 

a,(i) ~ (a,(0) + al(0)) = ^/2l^x,{Q), (10) 

where the hermitian operator xi = (a^ -|- a|)/\/2 is the 
field quadraure and (n^) = sf is the mode occupation. If 
we assume that highly populated modes have dominant 
contribution to the field operator, then 
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V2Mi»(0)?A,(r,t). (11) 



Within this approximation any product of operators 5 
depends only on Xi operators (in general it could also 
depend on pi = {a^ — aJ)/i\/2) and the n-th order corre- 
lation function, 

/9„(ri, . . . ,r„;t) = {S^ri,t) . . . (5t(r„, i)(5(r„, i) . ..S{ri,t)) 
= ({5t (n , i) . . . (r„ , t)S{rn , t) . . . <5(ri , t)}sym) , (12) 

where the symmetrization is performed with respect to 
Xi and Pi. Thus the evaluation of the quantum average 
in Eq. (|12|l is equivalent to the procedure in which we 
replace the operators i^'s and pi's with c- number vari- 
ables Xi and Pi and perform an integral of the product of 
5(r,t) = J2i \/'2{ni)xiipi{r, t) function with the Wigncr 
distribution of the vacuum state. 



p„(ri, . . . ,r„;t) = J d{x,} J d{pt}W{{x,}, {p,}) 
(5*(ri,i)...<5*(r„,t)<5(r„,t)...<5(ri,i) = (13) 
d{x,} [ d{p,}W{{x,},{p,}Mn,t)\\ . .\6{r,,,t)\'. 
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Here 

w{{x,}, {Pi\) - n ^ (-^? - p^) ■ (14) 

i 

The integrals over piS can be calculated explicitly and 
we obtain 

p„(ri, . . . ,r„;i) = j d{xi}Y\_^exp {~x1) x 

|<5(ri,i)r....|5(r„,<)|2 (15) 

We stress that any correlation function can be evaluated 
according to the procedure described above and conclude 
that 

P{{x,}) = \[l=e^v{-xl) (16) 

i ^ 

gives the probability distribution of quadratures' ampli- 
tudes XiS. 

Thereby we generated a scheme for calculating the sin- 
gle realization of the scattering experiment: we randomly 
choose the set of values Xi with probabilities P{{xi]) and 
construct the <5(r,t) function. To find density distribu- 
tion in a single realization we take I5(r, t)p. This result 
is analogous to that obtained in Here we see the 
analogy to the so-called Wigner stochastic method 0- 
In our method we randomly choose the amplitudes of 
the quadratures of the macroscopically populated modes 
and to calculate any moment of the field operator (5(r, t) 
we simply perform a statistical average over variables Xi . 
In contrary to the Wigner method, we do not introduce 
any external noise to the dynamics of the system. The 
detailed discussion will be given elsewhere. 

Now we apply the method derived above to plot the 
density of the hallo of scattered particles, in a single re- 
alization of the collision. Since the density of particles is 
measured after the free expansion of the atomic cloud, it 
represents the momentum distribution of atoms. Hence 
we show the results in the momentum space, |5(k, t)p, 
where k is a wave- vector of the scattered atoms. We 
present both the cut as well as the column density, i.e. 
/ dkz\S{k,t)\'^ , see Fig. ^ AH the data corresponds to 
the time t, long after the collision was completed. 

As shown in the angular size of the coherence of 
scattered atoms is proportional to 1//3. We expect that 
by decreasing /3, and thus increasing the range of coher- 
ence, at some stage we should observe clear anisotropy 
(spikes) in the distribution of scattered atoms. The spikes 
would be an effect of the bosonic stimulation of scatter- 
ing into highly occupied modes. On the other hand, as (3 
increases, the range of coherence decreases. Then, scat- 
tering of separate atoms is practically independent. In 
such a case we expect that all the momentum modes 
available to scattered atoms should be almost uniformly 
occupied. These predictions agree very well with the re- 
sults obtained in the numerical simulations. 



Figure^shows the cut for fc^ = as well as the column 
density of atoms scattered out from the condensate for 
three different values of /3. The ratio a/ (3, which is a 
measure of bosonic enhancement 0, is in all cases the 
same, and equals a//3 = 5. This value corresponds to the 
regime of strong bosonic enhancement, where condition 
of high mode occupation is fulfilled. For the fixed radius 
of the condensate one can vary the value of P by changing 
the mutual velocity of the colliding clouds. In Fig.^ case 
a) , the collision can be regarded as slow and we observe a 
clear speckle structure. For larger values of p (Figure^ 
and^) the distribution becomes more and more uniform 
and resembles an average over many realizations. 

Notice that the increase of (3 can be achieved also by 
changing the size of the clouds, which is proportional 
to number of atoms. That explains the experimental 
results obtained in They are consistent with our 

data. While their condensates where rather slow the size 
(3 X 10'' atoms) was very big, giving f3 ~ 100. In this 
regime we expect that after coarse graining introduced 
by the CCD camera resolution, the momentum density of 
scattered atoms would be highly uniform. This is ex actly 
what we see in experimental data shown in Fig. 2b of |l2l |. 
If one would intend to see the speckles, he should reduce 
the velocity of colliding condensates and/or reduce their 
size by decreasing the number of atoms. Alternatively 
one can use tighter traps. According to our estimates it 
should be feasible to reduce the value of f3 to reach the 
/? ~ 20 region. 

Due to the large number of atoms, even the data ob- 
tained in a single experiment has some attributes of 
the statistical average 0. It is illustrated in the last 
column of Fig^ In this case we fixed the angle be- 
tween ki = [(t)',9',Q] and k2 = [</>' + (j),e',Q], keep- 
ing their length equal to Q and calculated the average 
/ (5*(ki)(5(k2) dfi' for three different realizations (dashed 
lines) and compared it with quantum mechanical aver- 
age (i5^(ki)(5(k2)) (solid line). The agreement is striking 
and getting better with increasing value of (3 (when the 
distribution becomes more uniform). 

In summary we show that single realization of the ex- 
periment of two colliding condensates contains an inter- 
esting information and may differ substantially from the 
picture obtained after averaging over many realizations. 
Especially for slow collisions we predict the lumpy struc- 
ture in the momentum distribution of scattered atoms, 
with clear spikes even in the column density. The spikes 
found in this Letter are analogous to the speckles in the 
light beam produced by a multimode laser |l3l |. Even 
closer optical analogy can be found in recent studies of 
the parametric down conversion |l4l |. 
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2P03B4325, 1P03B14729, 1P03B14629, and the Polish 
Ministry of Scientific Research and Information Technol- 
ogy under Grant No. PBZ-MIN-008/P03/2003 (M.T.). 
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FIG. 1: (Color online): Momentum distribution and normalized density matrix of scattered atoms for three different sets of 
parameters (/?, a); a) (20,100), b) (40,200), c) (60,300). Left panel corresponds to a cut with fc^ = 0, center to a column density. 
The right panel shows the comparison of the normalized density matrix averaged over many realizations (solid line) with the 
three different cases (dashed lines) of the density matrix calculated for a single realization averaged over the full solid angle. 
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